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Abstract 

The equation of gluodynamics with modified integral kernel is solved up to 2NLO (next- 
to-next-to-leading order) and 3NLO (next-to-next-to-next-to-leading order) in perturbative QCD. 
A relationship between the nuclear factor Ns and the ratio of the cumulant and factorial moments 
of the multiplicity distribution is examined. A comparison of theoretical results and experimental 
data in p-Pb and Pb-Pb collisions at 200A GeV and 546A GeV is carried out. 

Introduction 

A process of the multiple particle production is of the greatest importance in high-energy 
particle interactions. To describe it, a system of coupled QCD equations has been elaborated [1- 
3]. In case of the running coupling constant, these equations can be solved only within the 
perturbative QCD [4, 5]. Usually one starts with a more simple case of gluodynamics, i. e. 
setting aside the quark interactions. It gives a chance to draw some conclusions about behavior 
of the moments of multiplicity distributions in gluon jets. Taking quarks into account does not 
significantly affect the results and leaves their qualitative features unaltered [6, 7]. Therefore the 
case of gluodynamics is considered in the paper. 

I, 2NLO and 3NLO perturbative solutions of the gluodynamics equation 

The following equation describes gluon interactions without quarks: 

G' {y) = \ dxK {x)yl{y + In x + ln(l - x))[G{y + In x)G(y + ln(l - x)) - G(y)] . 



cc 

Here G(y,z) = ^(l + z")P^- generating fimction of the multiplicity distribution, Pn is 

n=0 

parton multiplicity distribution in gluon jet, z - auxiliary variable; G'(y) = dG/dy,y = in(p6/Qo) 
= \n(2Q/Qo) - jet evolution variable, which determines the energy scale, p is initial parton 
momentum, ^ - the angle of divergence of jet (assumed to be small and constant), Q-the jet virtuality 



squared four-dimensional momentum, Qo = const; /oiy) ' (^siy) running coupling 

1 • 1 • • ^ ^ 2;r 51 ln(2v)^ -3. . , . . 

constant, wnicn is wntten as (y) = (1 ) + 0{y ) m two-loop approximation. 
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The moments of multiplicity distribution are defined as 

00 

F, = ^ = ^_^!^| „ - factorial, andi^, = ^ dHnGiz) ^ _ 

cumulant moment of q rank; (n(y)) - average parton multiplicity of jet with given j;. 

The kernel of the equation is ^(jc) = — -(1-jc)[2-jc(1-jc)]. The nuclear QCD factor Ns 

X 

N 

usually is inserted into the kernel in the following way [8, 9]: K(x) = — - - (1 - x)[2 - x(l - x)] . It 

X 

accentuates the change of ratio between the soft and hard parts of parton spectrum in the nuclear 
media. The singular part 1/x describes low-energy particles, and when the Ns > 1 emerges in 
numerator, the soft part of spectrum is spreading, i. e. the number of soft emitted gluons 
increases. It simulates the softening of the spectrum of produced particles due to their 
rescattering in the nuclear media. In our case it is convenient to choose the kernel as 

K(x) = — — ^— (l-x)[2-x(l-x)] . It doesn't change the physical interpretation of the 

X Ns 

phenomenon, but allows us to avoid rising of the coefficient ^Jn^ by the leading term of 
perturbative series. 

The dependence of mean multiplicity on energy is defined by QCD anomalous dimension 

as = • The lower limit of integration is constant and influences just the 

normalization factor, which is free adjustable parameter. That is why the renormalization of 
K(x) doesn't affect the result. For the QCD anomalous dimension y we use approximate 
expression / = 7o(l"~^i7o "^2/0 "^3/0 "•••) ~ high-energy (i. e. for large j;) perturbative series. 

To solve the equation in 2NL0 approximation, one should expand each function with 

2 

shifted argument in Taylor series and keep the terms up to -order in the right side of 
equation: 

^[In G(j;)]r = G(y)-l-2h,G'(y) + h,G''(y) + 4Bhy,[G(y)-l] . 

To 



From the last equation we get coefficients a i and a2, comparing the terms by the z by first 
power: 0^ = }\-^, = (^-h^Xa^ + lB)-^. Further, examining the z^-terms with arbitrary 

q, the explicit expression for Hg can be obtained: 

H = , l-^Vn + (M' + 2Vig + 4g/;iVo' f^nj^tion hag one minimum, located at 

" q' - {2a,q' + Bq)y, + [(af - 2a2)q^ + 2Ba,q]yl 

1 a 

^min ~ + ~ ' asymptotically tends to constant /jj/q for large q. When Ns increases, the 

first root of Hq moves to the left (in the direction of decreasing q), while the minimum position 
moves to the right. Fig. 1 shows the shift of the first root oiHq. 

In the 3NLO approximation (keeping the -order terms) we come to the following 
equation: 

' - [In G{y)-\r = G(y) - 1 - 2h,G'(y) + hfi\y) + 4Bhy,[G(y) + G'^y) + 



n 

+ h,2[G"iy)ilnG{y)y + G'(j)(lnG(j))"] - 4Bhy,G'iy), 

413 395 11;7"^ 

where = 2^(3) , = ^(3) + (^ - Riemann zeta funcfion) [11]. 
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This equation allows us to find the as coefficient 

3 1 5 1 17 

(flj = flj (flj - /?j + — 5) - — 55j + {a^ + —B)h2 - —h^, B^= — ) and gives a recurrence relation for 

calculation high-order Hg values via lower Hq 's using the well-known formula for cumulant and 
factorial moments = liC^^F^ . The graphs of Hq, plotted for Ns = 1,0.. 1,08 (fig. 2) and 

k=0 

Ns = 1,0.. 1,04 (fig. 3), demonstrate that when Ns is growing, the first root of Hq now moves to the 
right and amplitude of Hq oscillations decreases. For Ns > 1,2 disappearance of the oscillations 
may be observed, and behavior of Hq becomes the same as for 2NLO: one minimum and 
presence of asymptote. 



Comparison with simulated data 

In order to ascertain which approximation better corresponds with experimental data, the 
simulated multiplicity distributions for nuclear collisions of p-Pb and Pb-Pb at the 200A GeV 
and 546A GeV energy were examined. Using them, the Hq values were calculated and their 
graphs were plotted (fig. 4, fig. 5). One can see that the behavior of Hq is more precisely fitted by 



3NLO approximation, since the first root of Hq moves to the left when Ns is growing, and 
amplitude of oscillations decreases. Comparison of fig. 3 and fig. 4 permits to estimate the value 
of iV, for Pb-Pb collision: Ns ~ 1 ,0 for 200A GeV and Ns ~ 1 ,04 for 546A GeV. 
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Fig. 1. Shift of the roots of Hq according to 
increasing Ns in the 2NL0 approximation. 
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Fig. 2. Shift of the roots of Hq according to 
increasing Ns in the 3NL0 approximation. 
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Fig. 3. Behavior of Hq when Ns is growing 
in the 3NL0 approximation. 



Fig. 4. Hq for the collisions of Pb-Pb nuclei 
with energies 200A GeV and 546A GeV. 
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Fig. 5. Hq for the collisions of p-Pb nuclei with energies 200 A GeV and 546 A GeV. 



